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ABSTRACT 



We estimate the sum of the TBB couplings using QCD Spectral Sum Rules (QSSR). 
Our result implies the phenomenological bound £,'(vv' = 1) > —1.04 for the slope of the 
Isgur-Wise function. An analytic estimate of the (physical) slope to two loops within 
QSSR leads to the accurate value £,'{vv' — 1) ~ —(1.00 ± 0.02) due to the (almost) 
complete cancellations between the perturbative and non-perturbative corrections at the 
stability points. Then, we deduce, from the present data, the improved estimate \V c f,\ ~ 
(1.48ps/r B ) 1/2 (37.3 ± 1.2 ± 1.4) x 10~ 3 where the first error comes from the data analysis 
and the second one from the different model parametrizations of the Isgur-Wise function. 
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1 Introduction 



With the advent of the heavy quark symmetry there has been considerable interest 
and progress in the understanding of the semileptonic form factors of the transition of 
a heavy quark into another heavy quark, as in this infinite mass limit all semileptonic 
form factors reduce to the single Isgur-Wise function @. The question is whether this 
result in the heavy quark infinite mass limit can be applied to the physical B and D 
mesons. Progress has been made for including the 1/Mq corrections. At the non-recoil 
point, the 1/Mq terms cancel in the form factors ||, while for the leptonic decay constant, 
one expects a large 1/Mq correction as indicated by lattice j| and QCD Spectral Sum 
Rules (QSSR)P estimates. QCD vertex sum rules have also been applied for evaluating 
the semi-leptonic form factors at q 2 = @-[]8| and the Isgur-Wise function at (/mazllH 
iOj . Although the QSSR results are quite impressive and the authors ||-[]10| expect that 



working with the ratio of the vertex over the two-point pseudoscalar sum rules eliminates 
different systematic uncertainties, one should not forget that, in this derivation, certain 
assumptions on the choice of the QSSR scales and QCD continuum thresholds, which 
introduce uncertainties, have to be done. Most of these choices are based on the model 
discussed in Ref. |TT|. More rigorous and less model dependent is the bound obtained 



recently by ||12|| , where one exploits the analyticity and the unitarity of the 6-number form 
factor in the limit M c = Mb (see also [|H|,[|I4|]). The bound on the slope of the form factor 
in [12] depends on the residues and on the positions of T poles. The most rigorous bound 
comes from the normalization condition F(0) = 1 of the 6-number form factor leading to 
F'(v.v' — 1) > —6, while the inclusion of TBB couplings, assumed to be of order one 



each, leads to phenomenological bound [[12]j : 

F'(v.v' = 1) > -1.5, (1) 
which is 4 times much stronger than the previous one. 

The purpose of this note is to present an estimate of the TBB couplings which play 
an important role in the phenomenological derivation of the previous bounds ||12j| , |fL3| . 
We shall also estimate directly this slope from analytic expressions to two-loops in the 
QSSR approach. Then, we shall deduce an estimate of the CKM mixing angle V^,. 

As the BB pairs are below the T(1S,2S,3S) states, there are no available data for 
estimating such couplings. The only experimental available information concerns the 
coupling of the T(4S) which lies above the BB threshold. From the data of its leptonic 



and total widths, one can deduce [0 : 



where 



and 



\ VT4 \ = ^M< 0.75 ±0.15, (2) 
27t 4 



1 2 /l\ 2 M T4 



r(T(4^e+0 = -™^~j (3) 



r (t(4S) 



BB)=M rA ^)U-^'\ (4, 



with the normalization: 

<0|-&7' i 6|T i >=e^^-, (5) 
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where we have included in the definition of the current the negative sign related to the 
6-quark charge. In principle, the couplings of the three lightest states are unconstrained. 
However, using the location of the three T poles and the fact that the b-number form 
factor of the B-meson is 1 at q 2 = 0, Ref . | T!| obtains the following bounds: 



\VtA < 9.9 x 10 3 , |r/ Ta | < 17.1 x 10 3 , |r/ T3 | < 8.1 x 10 3 , (6) 

which, as already pointed by the authors, are quite weak compared to the one for 774 
obtained previously from the data. 



2 QSSR estimate of the TBB coupling 

In so doing, we consider the vertex correlator : 

V^q 2 ,p 2 ,p' 2 ) = {if [ d A x dV ( " ,?/) < 0|TJ 5 (a;)J'*(0)4(y)|0 > 



= (p + pTV + (q 2 ,p 2 ,p' 2 ), (7) 

built with the local quark currents: 

J 5 (x) = (m d + M b ) : d(t l5 )b : , J»= -: bfb : . (8) 
V + obeys the double dispersion relation : 

V + (q 2 y,p' 2 ) = -- 2 Hr^r T^lmV + (q 2 ,t,t') + ..., (9) 

7T Z JM% t — p z J Ml V — p u 

where ... means polynomial subtraction constants and q 2 = (p — p') 2 < 0. For the 
estimate of the VTi couplings, we evaluate the correlator at q 2 = similarly to the 
estimate of the form factors of the B semi-leptonic decays using vertex sum rules. We shall 
limit ourselves to the lowest order contribution in a s but include the non-perturbative 
condensate contributions in the OPE. This approximation has also given a quite good 
prediction in different estimates of the B and D decay form factors. So we, a priori, expect 
that a similar feature will hold in our analysis. The QCD expression of the three-point 
function has been evaluated in the literature ||, including higher dimension condensates. 
Our case corresponds to putting M b = M c in this paper. Therefore, the perturbative 
contribution reads, to leading order in a s : 

lmVr\q 2 , t, * -\q 2 - W - ^ (10) 

which shows that at q 2 = 0, the lowest order perturbative contribution to the spectral 
function vanishes. We have also checked this result from a direct evaluation of the triangle 
perturbative diagram. Therefore, in this case, the leading contribution comes from the 
light quark condensate and reads : 

n ( o,^) = ^<a> 1 (id 
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which is not the case of the other B decay form factors studied previously at q 2 = 0. 
So, from this particular feature, we (intuitively) expect that V + (0) is much smaller than 
the previous form factors. We parametrize the phenomenological side of the vertex by 
introducing the B and T couplings via: 



< Q| J 5 \B >= V2f B M- 



2 

Bi 



Ml 

< 0|J^|T, >= e^— ^ (12) 

and we insert the intermediate states in (7). Using the definition of the TBB coupling in 
(2), we have at q 2 = 0: 

(13) 

If one uses the quark hadron (semi) local duality picture by simply equating the phe- 
nomenological and QCD sides of the vertex, one obtains the sum rule: 

y. T ~r^ 3 ^^ d4) 



Using for f B the local duality constraint [15 



2f B Ms* 1 (#)'(£.)', (15) 



7T 2 V c ) \M B 

where ~ 1.3 GeV |5[] is the B continuum energy, one obtains : 

E Vr, * j^L < dd >^ -0.07, (16) 

if one uses < dd > (E^) ~ -(250MeV) 3 . The previous constraint indicates that r)y i 
remains constant for M^ — >• oo, in agreement with the expectation from the "naive Mf, 
counting rule": 

g TBB ^Ml 12 , 7 t^M 6 1/2 . (17) 

This feature increases our confidence on the physical meaning of the constraint in (14). 
Moreover, the previous constraint also indicates that the sum of couplings in (16) is almost 
independent of the b-quark mass as we shall check later on in the case of the J/0- 

One can improve further the previous constraint by working with the Laplace (Borel) 
sum rules and by including the contribution of the quark-gluon mixed condensate g < 
da^^-G^d >= Mq < dd >. In this way, the relativistic Laplace operator leads to the 
change : 

' re'™? , „ * ^r'e- T ' M K (18) 



p 2 -M 2 ' p' 2 -M 2 

where r (resp. r') are the sum rule variables associated top 2 (resp. p' 2 ). For convenience 
and because of the symmetry of the vertex, we shall take the natural choice r = r' in our 
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analysis. However, this choice does not have any noticeable effect in our conclusion. We 
use the following values of the QCD parameters pfl , [Pf] : 

f B ~ (1.59 ±0.26)/, 
M b (p 2 = Ml) ~ (4.59 ± 0.05) GeV 

< dd >= -((189 ±7) MeV) 3 (- log(r- 1/2 A))~ 2//31 

Ml ~ (0.80 ± 0.10) GeV 2 , (19) 

where -f3 x = |(11-^) for SU(n) f and A = 260±50 MeV. We give the result for J2r] Ti in 
Fig. 1 versus the sum rule variable r. Due to the absence of the perturbative contribution 
in our leading order analysis, the result is insensitive to the continuum contribution. The 
stability in r is reached at 0.2 GeV -2 , a value that is quite similar to the one where fs is 
optimal ||. At the stability point, one obtains: 

5> Ti ~ -0.224, (20) 

i 

where the error is negligible (0.004), if one uses the correlated values of js and Mb. 
We estimate the maximal error by taking the uncorrelated values of the previous two 
parameters. Then, we obtain, at the order we are working, the estimate : 

Vr t - -(0.224 ± 0.064). (21) 

i 

An improvement of this result needs an evaluation of the radiative correction to the 
vertex function. However, we do not expect that the higher order terms will modify the 
present leading order results by more than a factor 2-3, if no marginal terms break the 
conventional OPE or/and if there are no anomalous couplings that drastically modify the 
parametrization of the spectral function in (13). 

The result from the simple local duality relation in (16) corresponds to the case where 
r — > 0. The difference between the local duality and Laplace sum rules results can indicate 
the possible large role of the continuum for r — ► which is negligible in the Laplace sum 
rule analysis. 

We test the quality of our estimate by applying the method in the J /ip channel where the 
couplings of the 3S and 4S states are also bounded experimentally to be : 

Tfa < 0.47 , V ^ < 0.17. (22) 



As we have discussed in the case of the decay constant [0 , the simple duality constraint 
in (15) reproduces correctly the ratio /b//d, though its prediction for the absolute value 
is inaccurate. Then, we also expect that the ratio of the coupling ^ is well reproduced 
by (16). It gives : 

£^-(§|) £^ 0.40, (23) 
where we have used E® ~ 1.08 GeV p|,[15|. A direct Laplace (Borel) sum rule estimate 



anologous to the one used to get (21) gives an optimization scale r ~ 0.8 GeV 2 , very 
similar to the one for f D §. Using f D ~ (1.31 ± 0.12)/, and M c (p 2 = M 2 ) ~ (1.45 ± 
0.05)GeV 0, the sum rule gives the estimate : 

5>fc -0.34 ±0.02. (24) 



4 



This result compares quite well with the previous leading order estimate from local duality 
sum rule in (23) and with the experimental bounds in (22). It also indicates that the sum 
of couplings is almost independent of the heavy quark mass, which is manifest in the local 
duality constraints. This test increases further our confidence on the numbers obtained 
in (21). 

Moreover, our result applies if the quark structure of the vertex remains valid for this 
particular process but cannot be used if the BB, B*B* pairs are molecules formed by Van 
Der Vaals like forces. In this case a vertex sum rule approach with four-quark currents 



similar to the one done for the KK molecule [17| becomes more adequate 



3 Bound on the slope of the Isgur-Wise function 

We have estimated the sum of the TBB couplings using a conventional vertex sum rules 
analysis within the quark structure of the BB meson pairs. Despite the leading order 
approximation that we have used for deriving the values of the couplings, we expect that 
our results are valid within a factor of 2-3, which is a conservative estimate of the radiative 
corrections not included here. A comparison of the result with the experimental bound 
in (2) indicates that the bound for the T(4S) coupling is satisfied by the sum of the 
couplings of the different T states. In order to see the effects of our results for the bound 



on the slope of the Isgur-Wise function as derived in [12], we shall consider the following 
different scenarios on the strength of each coupling given the constraint in (19) for the 
sum: The first scenario, which seems to be the most plausible phenomenologically, due to 
the experimental suppression of the electronic widths of the 2S and 3S states, is the one 
where the coupling of the IS state is much larger than the previous ones (vector meson 
dominance) (first row in the tables). The second and third scenarios (second and third 
rows) are the ones where the absolute values of the couplings are equal, with also the 
possibility to have a cancellation between the couplings of the 2S and 3S. The fourth 
possibility (fourth row) is the one where the coupling of the 3S almost saturates the 
experimental bound for the 4S given in (1). The one of the 2S is assumed to be about 
the meson mass squared ratio (M3/M2) 2 times the one of the 2S. 

The upper and lower bounds on the slope of the 6-number form factor are given in these 
tables 0]. One can notice that increasing the sum of couplings by a factor 3 only affects 



the bounds by 15%. From the tables, one can deduce the conservative phenomenological 
bounds for the slope of the 6-number form factor: 

- (0.88 ~ 1.34) < F'(vv' = 1) < (0.08 ~ 0.52). (25) 

The lower bound is comparable with and even slightly stronger than the conservative 
bound in (1) given by Ref. [0, while the upper bound is weaker than the Bjorken bound 
of -1/4 []B|. Using its relation with the slope of the IW function £' [BCT]: 

e'(l)^F'(l)-^loga s (M 6 ), (26) 
75 

we can deduce the conservative lower bound: 

g(vv> = 1) > -1.04. (27) 
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to 


v / low&v 


F'(i) 


-0.224 


0.000 


0.000 


-0.967 


0.447 


-0.224 


-0.224 


0.224 


-0.882 


0.516 


-0.080 


-0.080 


-0.080 


-0.950 


0.465 


-1.760 


+0.890 


0.720 


-1.183 


0.235 



Table 1: Upper and lower bounds for the slope of the 6-number form factor for various 
phenomenological values of the couplings r/j with J^Vi = —0.224. 





m 




(1) hower 


F'(l) 

- 1 V / upper 


-0.672 


0.000 


0.000 


-1.118 


0.306 


-0.672 


-0.672 


0.672 


-1.175 


0.250 


-0.224 


-0.224 


-0.224 


-1.048 


0.373 


-2.170 


0.780 


0.720 


-1.344 


0.081 



Table 2: The same as in Table 1 but J^Vi = ~ 3 x 0.224. 



The previous results question the accuracy of the experimental domain [21|: 



- 2.3 < = 1) < -1.17, (28) 

obtained after extrapolating the data until the non- recoil point. Our results also indicate 
that the smallness of the sum of the TBB couplings, and presumably, of each individual 
coupling derived in this paper, raises again some doubts on the accuracy of existing models 
and methods used for determining the mixing angle V c f,. 



4 QSSR estimate of the slope of the IW function 

In view of the former result, let us estimate £'(1) analytically, using QSSR in the Heavy 
Quark Effective Theory (HQET). In so doing, we work as |fjfl-||10|, with the ratio of vertex 
over two-point function sum rules to two loops within the continuum model of |TTJ . Using 
the expressions in [I 



10], we can deduce the compact analytical expression of the physical 



Isgur-Wise function: 



<qq> T 



£phys(y = W') 



i + y, 



i + -m) 

7T / 




1 + v 2 ' 38 ) ( 1 ~ (its 



a 



—g(y) 



71 



+ < a s G 2 > t 4 



3/, 



+ < qq > t — — 




(29) 



6 



with: 

16 , „ 49\ „ / 8 , „ 197\ 

600/ 

.1.6 , , 5G\ . . ■■> / S , , 112\ 

225/ 



f(y) = l{y)h{r) + (y - 1) log 2 - H) - (1 - y) 2 (A log 2 - 
= 7 (y) (^(-^ c r) - lE ) + (y - 1) log 2 - p) - (1 - yf (A l og2 - 



(30) 



where: 



4 / log(y + VF^l) ^ 



7(?/) = 3 I V VW^T ~ 1 I ' 7(1) = °' 7 ' (1) = 8/9, 



lE = 0.5772, Ei{-x) = - I ^e~ x (31) 



CO 2 . . / /"^ Co" 2 



and: 

4(r) = P dw T -^ sr log(o;r)e— / I J = / <l^——< \ ■ ( 32 ) 

We shall use in our numerical analysis the value: u c ~ (3.0 ± 0.5) GeV with a generous 
error compared to the true error of 0.1 [pH , ||15|| . We leave r as a free parameter, which we 
shall fix from a variational method. The slope £'(1) is the value of the first derivative of 
the IW function with respect to y at y — 1, which we can deduce analytically from the 
previous sum-rule expression of £ p h ys - Then, its expression reads: 



£'(y = !) = -! + 6^ + 6 NP , (33) 



where: 



a s . . / ,. , 16 , 49 
<W * -^(r)^7'(l) + ylog2-- 

87r 2 o / a, 5 9 9 

8 NP ~ -— r 3 < qq > 1 + 1.05^ - -M V 

+ ^- 4l #^. (34) 
3i 967r 

We optimize this previous sum rule for £' using a variational method. The stability of the 
result is reached for r _1 ~ 1.7 GeV (which is about the characteristic scale of the b into c 
transition), while the result is very stable (less than 2% change) for a larger range of r _1 
between 1.25 and 2.5 GeV. An analogous stability is also obtained by moving uj c in the 
range between 2.7 GeV (starting of r stability) to 3.5 GeV. We also notice that there is 
an almost complete cancellation between the perturbative radiative and non-perturbative 
corrections where each strength does not exceed 4%. This feature leads to the accurate 
estimate: 

£'(1) ~ -(1.00 ±0.02). (35) 

We can multiply the previous error by a factor two in order to add a conservative system- 
atic error inherent in the method and in the continuum model. 

This value differs from the previous result in |J deduced from a numerical polyno- 
mial two-parameter fit of the Isgur-Wise function, where the errors given there are quite 
doubtful. Our result satisfies the bound derived previously and the one of [12]. 
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5 Improved estimate of V c ^ 



As an application of our result, let us estimate the mixing angle V cb by using the ARGUS 
pT| and CLEO [^] data on £(y) x V cb . We shall use in our numerical analysis the following 



parametrizations: 

e(y)~i+e'(y-i) 

£(<,)~exp{£'(y-l)} 
ay)^{(i + y)/2} 2e 

^)^(^ I )e XP {(^ + l)^|}. (36) 

The two former are in line of the Taylor expansion used for y around 1. The third is the 
pole parametrization and the fourth is based on overlap integrals of meson wave functions 
in a harmonic oscillator model. We shall normalize our result with the world average 
t b ~ (1.48 ± 0.10) ps of the .B-lifetime given in |2"3]] : 



^(i^W (37) 



tb 

Given the previous value of the slope in (35), we determine V cb from each data point and 
then, we make a weighted average of the different results. From B° — > D*~lv, we obtain 
in units of 10 3 : 

\V cb \ ~ 36.9 ±3.3 ARGUS 91 

~ 36.0 ±1.5 CLEO 93, (38) 

while from B — > D*°W, we obtain in units of 10 3 : 

\V cb \ ~ 40.6 ± 2.3 ARGUS 92. (39) 

Our results are the average of the ones from the previous alternative parametrizations 
of £ in (36). The last two parametrizations give almost the same results of V cb . The 
errors given there are the largest ones from each parametrization and are only due to the 
data. The choice of the parametrizations induces an extra error of 1.4. Theoretical errors 
induced by the ones of the slope are negligible. We take the average of previous results 
in (38) and (39). Then, we obtain the final best estimate: 

\V cb \ ~ (37.3 ± 1.2 ± 1.4) x 10" 3 , (40) 

where we the first error comes from the data, while the second one is due to the different 
choices of parametrizations used in the literature. We consider this result as a noticeable 
improvement over the existing estimate of V cb (see e.g. |23|], [p4j| ) thanks to a better 
control of the value of the slope from the QSSR estimate in (35). This result also agrees 
with the lesser accurate estimate (after rescaling the lifetime used in |7]]) from the B into 
D,D* semi-leptonic decays within a QSSR estimate of the form factors at zero momentum, 
where finite corrections due to the c and b quark masses have been taken into account. 

Moreover, a model-independent result from the phenomenological de Rafael- Taron-like 
bound in (27) gives: 

\V cb \ < 38.9 x 10~ 3 , (41) 
which is enough strong for eliminating some results given in the literature. 
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6 Conclusions 



We have estimated in Eq. (21) the sum of the TBB couplings using vertex sum rules. 
Using this information into the analysis of [121, we have derived in Eq. (27) a phenomeno- 
logical bound on the slope of the IW function. Finally, we have reestimated in Eq. (35) 
this slope analytically from QSSR. From the previous results, we have deduced the value 
of V c b in (40) and the bound in (41). The accuracy of this value of is mainly do to 
the good control of the slope both from the sum rules and from the de Rafael- Taron-like 
bound. 
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Figure captions 



r-dependence of the TBB coupling in the Laplace sum rule analysis. 
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